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1 Introduction 

■ Much of our knowledge about QCD and the structure of hadrons (mainly nucleons) 
' has been gained from Deep Inelastic Scattering (DIS) experiments such as eA^ eX 
or lyN jjL^ X. The (inclusive) cross sections are determined by structure functions Fi 
C\| and F2 when summing over beam and target polarisations (and an additional ^3 when 
^""^ using neutrino beams), and gi, 92 when both the beam and target are suitably polarised. 
Structure functions are functions of the Bjorken variable x and Q^, the large space-like 
momentum transfer from the lepton. (Another class of structure functions - the transver- 
sity hi - can be measured in Drell-Yan processes or certain types of semi-inclusive pro- 
cesses.) 

As a direct theoretical computation of structure functions does not seem to be pos- 
(— I sible, we must turn to the Wilson Operator Product Expansion (OPE) which relates mo- 
• • ments of structure functions to (nucleon) matrix elements in a twist (ie operator [dimen- 
sion - spin]) or Taylor expansion in So first defining bilinear quark operators 



X 
J3 



H' 0r;pi -p,. ^^pPi-A'.j^) ...iD q, (1) 



where q is taken to be either a u or d quark and T is an arbitrary Dirac gamma matrix we 

4) 



have for the nucleon matrix elements, the Lorentz decompositions (s^ = —rnj^) 



^ n+1 



s{'^pMi...pM.}_tr 



(s>{^i - p'^s^'"' • • • p^" > - tl- 
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{mi _ r>'^e{A'n„'^2 . . . 



(2) 



where the symmetrisation/anti-symmetrisation operations on the operator indices also 
indicates that they are traceless (which gives them a definite spin). w„, a„, d„ and t„ can 
be related to moments of the structure functions. For example we have for w„ and F2 

dxx-'F,(x, Q') = l Eg^^ip^Q', gnvi^^'^if^) + 0{1/Q') , 

f=uA,g,... 

(3) 

and similar relations hold between gi and a„; 32 and a linear combination of a„ and d„; 
hi and Although the OPE gives t;„ from Fi (or F2) for n — 2,4,...; w„ from F3 
for n = 3, 5, . . .; a„ from 51 for n = 0, 2, . . .; a„, (i„ from 52 for n = 2, 4, . . ., other 
matrix elements can be determined form semi-exclusive experiments, for example ai by 
measuring tt* in the final state. 

While the Wilson coefficients, E^''^{1, g'''^{Q)) are known perturbatively (typically 
two to three loops) and determine how the moments change with scale, the 'inital con- 
dition' ie the matrix element is non-perturbative in nature. The only known way of 
determining them from QCD in a model independent way is via Lattice Gauge Theory 
(LGT). In this talk we review our status (QCDSF and UKQCD Collaborations) of some 
aspects of these determinations including some higher twist results. A more general re- 
view may be found in (Gockeler et al 2002a). We shall also restrict ourselves here to 
forward matrix elements (and so not consider the form factors and the more embracing 
Generalised Parton Densities or GPDs). Determining moments of structure functions is 
an active field of reseach at present, see for example (Dolgov et al 2002, Guagnelli et al 
2004, Ohta et al 2004). 



2 The Lattice Approach 

The lattice approach involves first Euclideanising the QCD action and then discretis- 
ing space-time with lattice spacing a. The path integral then becomes a very high di- 
mensional partition function, which is amenable to Monte Carlo methods of statistical 
physics. This allows ratios of three-point to two-point correlation functions to be de- 
fined, 

= {N.{t,p^O,ij)N,i,;P)) ^ ^^^^^^^o,\N,m , (4) 
{N{t;p)N{0;p)) ^ " 

where Na is some suitable nucleon wavefunction (with Dirac index a) such as 

X 

The proportionaHty holds for ^ r <C i < ^Nt for a lattice of size Ng x Nt- There 
are two basic types of diagrams to compute in eq. ©: the first is a quark insertion in one 
of the nucleon quark lines ('quark line connected'), while in the second type the operator 
interacts only via gluon exchange with the nucleon ('quark line disconnected'). Due to 
gluon UV fluctations these latter diagrams are numerically difficult to compute. However 
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by considering the Non-Singlet, NS, or Ou — Od operators, giving matrix elements such 
as Vn-NS = i^n"^ — v^'' in eq. ^ then the f ^ s and g (gluon) terms cancel. (For 
higher moments however, one might expect that sea effects anyway are less significant 
as the integral is more weighted to a; ~ 1.) Although LQCD is in prinicple an 'ab 
initio' calculation there are, of course, several caveats. First our lattice 'box' must be 
large enough to fit our correlation functions into. A continuum limit a ^ must be 
taken. A chiral extrapolation must be made from simulations often at the strange quark 
mass or larger down to the almost massless u/d quarks, or until we can match to Chiral 
Perturbation Theory, xPT (the problem there being that the radius of convergence of xPT 
is not known). Also to save CPU time, the fermion determinant in the action (representing 
n f quark flavours) is often disgarded - the 'quenched' approximation. Finally in addition 
to all the above problems the matrix element must be renormalised, in order to be able to 
compare with the phenomenological A/5 results. 

To attempt to address some of these issues we have generated data sets (Bakeyev et 
al 2004) 

1. 0(a)-improved Wilson fermions ('clover fermions') in the quenched approxima- 
tion at three couplings (3 = = 6.0, 6.2 and 6.4 (Gockeler et al 2004) corre- 
sponding to lattice spacings ~ 2.12, 2.91 and 3.85 GeV. (This checks lattice 
discretisation errors, which should be 0{a^).) The pseudoscalar mass, nips, lies 
between 580 Me V and 1200 Me V. 

2. Unquenched clover fermions at irips down to ^ 560 MeV in order to see if there 
are any discernable quenching effects. Various couplings are used, f3 — 5.20, 5.25, 
5.29 and 5.40 with lattice spacings ranging from a^^ ^ 1.61 GeV to 2.4 GeV. 

3. Wilson fermions at one fixed lattice spacing, a^^ ^ 2.12 GeV in the quenched 
approximation at pseudoscalar masses, nips, down to ^ 310 MeV, (nips /my ~ 
0.4) in order to try to match to chiral perturbation theory. This lattice fermion 
formulation has discretisation errors of 0{a). 

4. Overlap fermions, in the quenched approximation at one lattice spacing a^^ ^ 
2.09 GeV down to nips of about 440 MeV. These have a chiral symmetry even 
with finite lattice spacing and hence have better chiral properties than either Wilson 
or clover fermions (and also have discretisation errors of 0{a^)). 

Note that the physical pion mass is about ~ 140 MeV and we use the force scale 
ro = 0.5 fm = (394.6 MeV) to set the scale. These results cover various patches 
of {irips, a,nf) space. This is however not completely satisfactory. Overlap fermions, 
although the best formulation of lattice fermions known, are very expensive in CPU time, 
and are only just beginning to be investigated, eg (Galletly et al 2003, Ohta et al 2004), 
to which we refer the reader to for more details. 

The results obtained from eq. (0} are, of course, bare results and must be renor- 
malised. We shall not discuss this further here, just noting that many one-loop per- 
turbative results are known; but are generally not very satisfactory as lattice perturbative 
series do not appear to converge very fast. (The convergence can be helped using tadpole- 
improvement.) A preferred non-perturbative method is also available, (Martinelli et al 
1995) and the results presented here will have the Zs determined by this method. 
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All our results are for hadrons containing light (ie u/d) quarks. Reaching this limit 
is extremely costly in CPU time (and except for the overlap formulation, other problems 
connected with the non-chiral nature of the fermion formulation may arise). Much work 
has been done recently on chiral perturbation theory and it would be highly desirable to be 
in a region where these results can be matched to lattice results and then the limit nips 
TTiTT can be taken. Although one should take the continuum and chiral limits separately 
(and preferably in that order) we shall try here a variant procedure of using a simultaneous 
'plane' fit containing both limits. This is because at present the unquenched, data set 2 is 
less complete than the quenched data set 1 and this procedure at least allows for a direct 
comparison of results. (For set[2these different fit procedures lead to similar results.) 
Practically we might thus expect that for a quantity Q of interest 

Q = FQ{romps)+df{a/ror. (6) 

(romps) describes the (chiral) physics and the discretisation errors are 0{a'') where 
s = 1 for Wilson fermions and s — 2 for clover fermions. Naively one might expect a 
Taylor series expansion for F^ to be sufficient, ie 

F^{x)^Q{0)+c^x^ + ... , (7) 

where x = romps- Over the last few years expressions for F^ have been found 

fQ{x) = Q(0) (1 - cQx' \n{x/roA^)') +..., (8) 

showing the existence of a chiral logarithm ~ Inmg (including the quenched case). 
For Vn;NS^ dn-NS, tn;NS the Constant is known (and positive), see eg (Chen et al 
1997). One expects most effect of the chiral logarithm for ti-^^s and least for Qo-ns- The 
chiral scale, A^, is usually taken to be ^ 1 GeV. The range of validity of the expansion, 
eq. (|8} is not known; one might expect that for mps > A^, pion loops are suppressed, 
leading to a smooth variation in ie constituent quarks, while for mps < A^ non-linear 
behaviour would be seen. Thus building in some of the constituent or heavy quark mass 
expectations, an equation of the form 

has been proposed (Detmold et al 2001). 

Present (numerically) investigated matrix elements include V2 = (x) (also part of the 
momentum sum rule: '''''' + (x)''^^ — 1), = (x^), = (x"^), oq = 2Aq (also 

occurring in neutron decay, as well as the Bjorken sum rule, as Au — Ad — and 
connected with quark spin), ai = 2Aq'-^\ 02 — 2Aq^^\ to = 2Sq, ti — 2Sq^-^^ and ^2- 
We shall only discuss w„, n = 1, 2, 3, ao/2, to/2 and ^2 here. 

3 Results of Continuum/Chiral Extrapolations for some 
twist two operators 

We now show some results. We start by considering w,'^*^ (2GeV) for n ~ 2, 3, 4. In 
Fig. n we show V2.fi g from data set 1, together with a fit using eqs. (|6j and Q. We see 
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Figure 1. WjfjvsC^GeV) versus (rQUips)^ (upper plot) and versus (a/ro)'^ (lower plot) 
using data set^\ Filled circles, squares and diamonds represent the three lattice spacings 
corresponding to j3 — 6.0, 6.2, 6.4. The chiral limit {ronips)'^ — is shown as a short- 
dashed line, while the physical pion mass is denoted by the long-dashed line. Also shown 
as a dot-dashed line is the mass of a hypothetical ss meson calculated as ^ \/2mK- The 
MRST phenomenological value is denoted by a star. 
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that O(a^) discretisation errors are small and seem to be relatively benign. By this we 
mean that the only limiting factor with the extrapolation is the amount of data available. 
We shall (thus) in future assume that this limit is not a problem. This does not seem 
to be the case with the chiral extrapolation where the data seems to strongly favour a 
linear extrapolation rather than the xPT result in eq. (|8|l. The value found in the chiral 
limit is about 50% larger than the MRST phenomenological value, (Martin et al 2002). 
(Note however that there are exciting hints that overlap fermions may be closer to the 
phenomenological value, (Galletly et al 2003, Giirtler et al 2004).) 

The same situation persists for the higher moments v^'^ and v^'^ . In Figs.|2]and|3lwe 
show these moments and compare the results with the MRS phenomenological values. 
In all cases we find that the moments are too large in comparison with phenomenological 
result. It is not clear why this is so, again the quarks seem to be acting more like con- 
stituent quarks rather than current quarks. Possible causes are quenching and/or a chiral 
extrapolation from too heavy a quark mass. We first consider possible quenching effects. 
In Fig. H we consider ^2^^ s '^S^^^' this time for n f — 2 flavours using data set 2. 
No real difference is seen in comparison to the quenched case. Indeed for other matrix 
elements considered a similar situation prevails. 
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Figure 2. ^35^5(2 GeV) versus (rgmps)^. Same notation as in Fig.^ 
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Figure 3. vf.f^g{2 GeV) versus (romps)^. Same notation as in Fig.^ 



To try to examine the situation at smaller quark mass, we now turn to the data set 3. 
Most of the above results have a quark mass at the strange quark mass (or heavier). In 
this data set we have generated quenched Wilson data at one lattice spacing, at light pion 
masses down to 310 MeV. In Fig.|5lwe show v^'.^gi^ GeV). In comparison with the pre- 
vious pictures note that the x-scale only runs to (rgmps)^ ~ 3.0. Again (except possibly 
for the lightest pion mass) the data seems rather linear (and constant). Also shown is a 
forced fit from eq. (|9j, leaving and free but constrained to go through the MRST 
phenomenological value at rups = mjr. Ignoring the lightest quark mass point, this is 
just possible; however it is very unnatural giving, for example, A^ ~ 500 GeV which is 
a very low value. 

A similar situation holds for the axial ao;jvs/2 = gA and tensor charge to.Ns/2. In 
Figs.|6land0 Again the results in both cases seem very linear Indeed from eq. (|9} due 
to the negative sign, we must have the lattice data decreasing to the phenomenological 
value. This is certainly not the case here (although experimentally tQ.jys is not known, we 
expect a similar situation as for g^; indeed in the non-relativisitic limit to-Ns/'^ ~^ 9a)- 
Later work including the A as well as the N in chiral perturbation theory (Detmold et 
al 2002, Hemmert et al 2003) reduce the c|-* coefficient, but this still is a problem. We 
show only linear fits, giving for gA a value somewhat lower than the experimental one. 
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Figure 4. V2'pfg{2GeV) versus (rQirips)'^ for unquenched fennions using data set 2. 
P — 5.20 results are (filled) circles; 5.25 squares; 5.29, diamonds; 5.40 down triangle. 
Otherwise the same notation as in Fig. 
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Figure 5. 772^^5(2 GeV) versus {rompsY for Wilson fennions from data set 3. The fit' 
uses eq. Otherwise the same notation as in Fig. 



(There are two possible caveats: for clover fermions the continuum extrapolation may 
have significant O(a^) effects, in distinction to Vn-NS and also there may or may not be 
larger finite volume effects present both in the data and theoretically, see for example the 
discussion in (Cohen 2001).) But at present the same general picture emerges as for the 
unpolarised moments. 



4 Some higher twist operator results 
4.1 Twist three 

The prime example is given by d2, which can be determined from §2, the first moment of 
which is a linear combination of 02 and d2. The operators for the a„ moments have twist 
two, but dn corresponds to twist three and is thus of particular interest. A 'straightfor- 
ward' lattice computation, (Gockeler et al 1996), gave rather large values for dj^ (where 
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Figure 6. ao-Ns/'^ ~ 9a versus {rQirips)^ for Wilson fermions (data set 3), with a linear 
fit. Otherwise the same notation as in Fig. 
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Figure 7. io;JVs/2 versus {romps)'^ using data set 3. The same notation as in Fig.^ 



= Q(")2^^"^ + Q'^'^'l^d'^^). A recent experiment, (E155 Collaboration 1999), how- 
ever indicated that this term was very small, which would mean that 52 is almost com- 
pletely determined by gi (the Wandzura-Wilczek relation). This problem was traced in 
(Gockeler et al 2000a), to a mixing of the original operator with a lower-dimensional op- 
erator. This additional operator mixes cx 1/a and so its renormalisation constant must be 
determined non-perturbatively. In (Gockeler et al 2000a) this procedure was attempted, 
and led to results qualitatively consistent with the experimental values. Note that this 
is only a problem when using Wilson or clover like fermions, as we would expect the 

additional operator to appear like ~ m^gcr Dq and hence vanish in the chiral limit. Thus 
there should be no mixing if one uses overlap fermions. 



4.2 Twist four 

Potential higher twist effects are present in the moment of a structure function, see eq. (|3}- 
These 0(1/ Q^) terms are composed of dimension 6, four quark matrix elements. A gen- 
eral problem is the non-perturbative mixing of these operators with the previous dimen- 
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sion 4 operators. At present results are restricted to finding combinations of these higher 
twist operators which do not mix from flavour symmetry. For the nucleon the S'C/i?(3) 
flavour symmetry group must be considered, ie taking mass degenerate u, d and s quarks, 
(Gockeler et al 2002b) giving 

dxF,{x, Q')\fJrHL.nn = -0-0005(5) ^^^^^ + Oial) , (10) 

for quenched Wilson fermions (ie part of data set 3). To access this moment experimen- 
tally needs very exotic combinations of moments from the measurement of the p, n. A, 
S and S baryons and is not possible. Nevertheless this term is very small in comparison 
with the leading twist result, and might hint that higher twist contributions are small. 

5 Miscellaneous Pion and Lambda results 

Moments for the pion and rho structure functions were computed in (Best et al 1997), 
for unimproved Wilson fermions. Using the Schrodinger Functional method, was 
recently calculated for the pion, (Guagnelli et al 2004) for both unimproved and 0{a)- 
improved fermions. A higher twist matrix element for the pion has also been computed 
for quenched Wilson fermions (ie using part of data set 3) 

dxF2{x, Q')\j,=lhtmann = 1.67(64) ^^^^^ + 0{al) , (11) 

where the SUf{'2) flavour symmetry group gives the combination = i^^^ +F2 — 
2F2 . This is again a rather small number, so although a rather exotic combination of 
pion matrix elements it might indicate that higher twist terms are small. 

Finally there have been results for moments of A structure functions, (Gockeler et al 
2002b) again using Wilson fermions (ie part of data set 3). These are potentially useful 
results as one can compare with nucleon spin structure and check violation of SUf{3) 
symmetry. First indications are that there is no evidence of flavour symmetry breaking in 
the matrix elements ie that A and p are related by an SU{3) p flavour transformation. 

6 Conclusions 

Clearly the computation of many matrix elements giving low moments of structure func- 
tions is possible. We would like to emphasise that a successful computation is a fun- 
damental test of QCD - this is not a model computation. There are however many 
problems to overcome: finite volume effects, renormaUsation and mixing, continuum 
and chiral extrapolations and unquenching. At present although overall impressions are 
encouraging, still it is difficult to re-produce experimental/phenomenological results of 
(relatively) simple matrix elements, eg V2. But progress is being made by the various 
groups working in the field. For example in comparison to our previous results, (Gock- 
eler et al 1996) there are now non-perturbative Zs and considerations of both chiral and 
continuum extrapolations and some unquenched results are now available. While the 
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continuum extrapolation seems to be 'just' a matter of more data at smaller lattice spac- 
ing, the chiral extrapolation does seem to present a problem, with no sign of any chiral 
logarithms being seen as predicted by xPT. Clearly everything depends on the data and 
the quest for better results should continue. To leave the region where constituent quark 
masses give a reasonable description of the data (ie linearity) unfortunately requires pion 
masses rather close to the physical pion mass. In this region fermions with better chiral 
properties will probably be needed, such as overlap, which in turn will need much faster 
machines. 
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